Strong Magneto-elastic coupling in VOC1: Neutron and synchrotron powder x-ray 

diffraction study 
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We present a combined neutron and synchrotron powder diffraction study of the crystal and 
magnetic structure of VOC1. The occurrence of antiferromagnetic order in VOC1 is accompanied 
by a monoclinic distortion. The high sensitivity of the magnetic interaction parameters on V-O-V 
bond angles and V-V distances yields strong magnetoelastic coupling which implies the structural 
distortion below Tat as well as anomalous structural effects in the paramagnetic phase. 



I. INTRODUCTION 

The transition metal oxyhalides MOX with M = Sc, 
Ti, V, Cr, Fe and X = CI, Br have attracted strong in- 
terest in the past years initiated through the challenging 
low-dimensional phenomena observed for TiOGi*2i2iSiiii£. 
At Ti = 67 K, TiOCl shows a transition associated with a 
dimerization of the Ti-ions along the 6-direction (first or- 
der phase transition) which has been interpreted as a spin 
Peierls transition. Thus, TiOCl might be a rare case of 
an inorganic spin Peierls system similar to CuGeOuii^i 
The whole group of transition metal oxyhalides MOX 
provides an interesting class of - at room-temperature 
- iso-structural compounds in which several electronic 
configurations and their interactions with the crystal lat- 
tice can be studied systematically. However, the physical 
properties of the MOC1 compounds different from TiOCl 
have been studied in less detail so far— ^ i 10 ' 11 ? 12 . Like 
TiOCl, VOC1 crystallizes in an orthorhombic structure 
with space group Pmmn, with a single crystallographic 
V-sitei 3 -. The V 3+ -ions are coordinated by four oxy- 
gen atoms and two chlorine atoms forming strongly dis- 
torted VO4CI2 octahedrons; see Fig. [2j These are linked 
via corners and edges forming quasi-two-dimensional bi- 
layers in the afo-plane. Due to the CI- ions there is a 
clear separation of neighboring Ti-O-bilayers explaining 
the electronic and magnetic low-dimensional character. 
Due to electronic correlations all oxyhalides are Mott- 
insulators. In VOC1 the band gap is reported to be about 
1.7 eV— . Whereas a direct overlap of the ground state 
yz-orbitals (with reference to the cell axis b and c) along 
the crystallographic 6-axis leads to one dimensional mag- 
netic interactions in TiOCl^, VOC1 is believed to exhibit 
quasi-two-dimensional antiferromagnetic ordering below 
Tj\r (80 K) 9 . In an early powder neutron diffraction study 
Wiedenman et al. reported a collinear magnetic order 
with the wave vector (1/2 1/2 1/2) and the magnetic mo- 
ments pointing in et-direction. The magnetic susceptibil- 
ity was well described within a quadratic 2D-Heisenberg 
model. However, no signs of a structural transition in 
VOC1 have been observed at low temperatures^. 

In this work, we report x-ray and neutron powder 
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FIG. 1: (color online) Synchrotron radiation powder x-ray 
diffraction patterns, (a) Rietveld fit (gray) to the data at 
15 K (black) with an enlargement of the diffraction pattern 
around the (±2 1 2) profiles in the inset, (b) Temperature 
dependency of the (±212) and (±121) profiles. 



diffraction experiments on VOC1 which reveal a struc- 
tural phase transition accompanying the antiferromag- 
netic order in VOC1. We find a remarkably strong magne- 
toelastic coupling in VOC1 which seems to arise from the 
dependence of the magnetic interaction on soft lattice de- 
grees of freedom, V-O-V bond angles and V-V distances, 
in analogy to the spin-Peierls compound CuGeO^ i 15 ' 16 . 



II. EXPERIMENTAL 

VOC1 was prepared by the chemical transport tech- 
nique as described by Schafer et al»ii. Single crystals 
obtained via this route were thoroughly ground, com- 
pressed under hydrostatic pressure and reground after- 
wards. This procedure was repeated several times in 
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order to achieve a homogenous powder and hence min- 
imize preferred orientation effects. The phase purity of 
the sample was checked several times using powder x-ray 
diffraction. 

Synchrotron radiation powder x-ray diffraction mea- 
surements have been performed at beamline B2 at Ha- 
sylab/DESY in Hamburg, Germany^. The single phase 
VOCl-powder was finely ground and filled in a glass cap- 
illary with 0.3 mm in diameter. The incident X-ray wave- 
length was 0.7466 A and the reflections have been mea- 
sured using an image plate detector in Debye-Schcrrcr 
geometry^. 

Powder neutron diffraction measurements have been 
performed at the D1A, D2B and D20 diffractometer at 
the ILL in Grenoble, France (A = 1.91, 1.594 and 2.4233 
A, respectively). The samples measured at the D20 and 
D2B diffractometer contained some unidentified impu- 
rity phases. Therefore, only the single phase sample 
measured at the D1A diffractometer was used for a full 
nuclear-structure refinement. The high-intensity D20 
data have been used for the magnetic structure refine- 
ment, as there was no overlap of peaks from the impurity 
phases with the magnetic peaks due to the longer wave 
length. 
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FIG. 2: (color online) Nuclear and magnetic structure of 
VOC1; yellow: C1-, blue: 0-, magenta: V-ions; black: the 
structural unit cell. The directions of the magnetic moments 
at each V-site is indicated by magenta arrows. Green (red) 
dotted lines indicate the antiferromagnetically (ferromagneti- 
cally) coupled V-V pairs in the AFM1 structure. 



III. RESULTS AND DISCUSSION 

A. Monoclinic distortion 

The temperature dependent powder x-ray diffraction 
measurements using synchrotron radiation clearly reveal 
a structural phase transition below 80 K. A diffraction 
pattern at the lowest temperature, 15 K, is shown in Fig. 
[T] (a) . The open circles denote the measured data points 




\ 

V. 


(a) 


90.2 ^ 
90 1 °*> 






90.0 






(e) 






/ 


98.5 


0.83% 




98.0 y 


~\ 


(f) 


0.068 

to 

0.067 



100 200 300 100 200 300 

T (K) T (K) 



FIG. 3: Results of synchrotron radiation powder x-ray diffrac- 
tion measurements: (a-c) lattice parameter, (d) monoclinic 
angle, (e) unit cell volume and (f) orthorhombic splitting 
e = (a-b)/(a+b). (Lines are guide to the eyes.) 



and the line indicates the Rietveld fit of the structure 
model using the program FullProf/WinPLOTR 18 . The 
part of the diffraction pattern around the (±212) reflec- 
tions is shown in the inset. The clearly observable split- 
ting of the orthorhombic (212) reflection indicates the 
symmetry reduction from the orthorhombic (Pmmn) to 
a monoclinic symmetry. The data can be well described 
in space group P2/n (monoclinic setting c) and a mon- 
oclinic angle 7 of 90.2°. In Fig. [1] (b) the (±212) and 
(±121) reflections are plotted for different temperatures 
indicating the orthorhombic to monoclinic phase tran- 
sition below 80 K. The resulting lattice parameters are 
plotted as a function of temperature in Fig. [3] (a-d). A 
rise of the monoclinic angle 7 can be observed below 80 
K reaching its maximum value of about 90.2° around 40 
K. The monoclinic distortion and the structural phase 
transition apparently have escaped detection in the pre- 
vious neutron diffraction experiments due to their limited 
resolution 9 . The group-subgroup relationship between 
Pmmn and P2 /n as well as the small change in the unit 
cell volume suggest a second order phase transition. 

Powder neutron diffraction measurements have been 
performed at the ILL. In Fig. 0] (a-b) neutron diffrac- 
tion patterns of VOC1 are shown. Fig. 2] (a) shows the 
magnetic peaks measured at 2 K at the high-flux D20 
diffractometer with an incident wavelength of 2.4233 A. 
In Fig. 0] (b) an entire diffractogram measured at 10 K 
on the DlA diffractometer using a neutron wavelength of 
1.91 A and higher resolution is shown. The latter data 
was used for structural refinements. 

The appearance of magnetic peaks at (1/2 1/2 1/2) 
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and (1/2 1/2 3/2) with respect to the nuclear lattice 
clearly indicates the antiferromagnetic order occurring 
below T^r ^80K in perfect agreement with Ref.— . A 
fully three-dimensional (3D) antiferromagnetic structure 
could be fitted to the neutron data; see Fig. [4] (a). The 
half-indexed reflections indicate the magnetic propaga- 
tion vector k = (1/2 1/2 1/2) corresponding to a 2 x 2 x 2 
enlargement of the magnetic unit cell. 

Within the nuclear unit cell there are two magnetic 
sites, the V ions at (0.75,0.25,0.117) and at (0.25,0.75,- 
0.117) which belong to two different chains of edge- 
sharing VO4CI2 octahedrons along the 6-direction, la- 
belled 6-chains in the following. These 6-chains form 
the VO-planes parallel to the a, 6-planes and are situ- 
ated alternatingly above and below the z=0 level. The 
inspection of the magnetic paths in Fig. 2 (a) and (b) 
allows one to qualitatively understand the antiferromag- 
netic order. There is a magnetic interaction J a along a 
mediated through the common O between two V sites 
at the same z-level. This antiferromagnetic interaction, 
which is absent in TiOCl due to the single occupation of 
the d yz orbital, only acts on V-ions in next-nearest neigh- 
bor 6-chains yielding an antiferromagnetic arrangement 
between these next-nearest 6-chains and thus a doubling 
of the a parameter. There is also an interaction Jf, along 
6 mediated through the common octahedron edges which 
results in an antiferromagnetic order within the 6-chains 
and doubling along 6. With an additional weak anti- 
ferromagnetic coupling perpendicular to the planes one 
may thus explain the doubling of the magnetic lattice 
along the three orthorhomic directions. However, in the 
orthorhombic phase the interaction between the nearest- 
neighbor 6-chains is fully frustrated since the path par- 
allel to [110] and that parallel to [1-10] are equivalent by 
symmetry and thus Jdiag = ^diag- Since these two interac- 
tion parameters connect to spins of opposite sign in the 
nearest-neighbor 6-chain the coupling is fully frustrated 
in the orthorhombic phase. 

The two possible magnetic arrangements in the or- 
thorhombic phase, however, do not correspond to two 
distinct magnetic symmetries but only to two domains 
of equivalent magnetic symmetry. This conclusion is 
fully corroborated by the representation analysis of 
the orthorhombic structure which only yields one irre- 
ducible representation for the propagation vector k = 
(1/2 1/2 1/2) and magnetic moments aligned parallel to 
the a, 6-planes. The monoclinic distortion accompanying 
the magnetic order lifts the frustration of the nearest- 
neighbor 6-chain coupling, J^ iag 7^ J^iag > smce ^ renders 
the magnetic paths parallel to [110] and [1-10] inequiva- 
lent. As indicated in Fig. 2 (b), all the Jdi ag -paths paral- 
lel to [110], Jdiagi are equivalent in the monoclinic phase 
but they differ from the J<ii a£ ,-paths parallel to [1-10], 
J diag- Supposing that J> ag is larger than J< ag the mag- 
netic ordering shown in Fig. 2 (b) is stabilized. We may 
thus conclude that the monoclinic distortion accompa- 
nying the magnetic transition lifts the frustration of the 
magnetic coupling between nearest-neighbor 6-chains. 
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FIG. 4: (a) Powder neutron diffraction patterns measured at 
2 K (D20 diffractometer) ; circles: data points, lines: Rietveld 
fits to the data with the AFM1-AFM4 structure models. The 
area below the AFMl-fit (moment \\a) is filled (gray), (b) Fit 
of the D1A measurement at 10 K; black: data points, gray: 
Rietveld fit to the data, It. gray: difference I f, s -Ica!c; arrows 
indicate the magnetic peaks. 



B. Magnetic structure 

Let us now turn to the quantitative analysis of the 
magnetic structure in VOC1. The two strong magnetic 
Bragg peaks observed at (1/2 1/2 1/2) and (1/2 1/2 3/2) 
are best described by the model with ferromagnetic cou- 
pling between the two V-moments at (0.75,0.25,+z) and 
(0.25,0.75,-z) pointing in a-directipon; see Fig. [2] Other 
directions of the magnetic moment like the [110]-direction 
induce discrepancies in the measured intensities of the 
(±1±11) and (±1±13) reflections (in notation of the 
magnetic unit cell). 

Due to representation theory, there are four irreducible 
representations for space group Pl\2/n with a magnetic 
ion at (-0.25 0.25 zq) (zq « 0.117) and with propagation 
vector k = (1/2 1/2 1/2). The basis vectors for these 
four possibilities derived from the program BasiRep 18 
are shown in Fig. 5. These four possibilities are re- 
ferred to as the antiferromagnetic structures AFM1 to 
AFM4. The representation analysis separates the ar- 
rangements with the moment parallel to c (AFM3 and 
AFM4) from those with in-plane moments (AFM1 and 
AFM2). For these two cases, there is the possibility of 
ferromagnetic or antiferromagnetic coupling of the two 
sites at (0.75,0.25,+z) and (0.25,0.75,-z). The two mag- 
netic structures AFM1 and AFM2 (or AFM3 and AFM4) 
can not be distinguished if the crystal structure is or- 
thorhombic (Pmmn). However, the structural phase 
transition Pmmn — ► P2/n observed in our work in- 
duces a splitting of the (±1 +1 1) and the (±1 ±1 3) 
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reflections (referring to the magnetic cell). The small 
splitting of these reflections allows us to distinguish be- 
tween the two magnetic structures. In Fig. |4] (a) the 
Rictvcld fits to the measured data at 2 K are shown 
for the four magnetic models. In the Rietveld fit of 
the AFM2 structure, the calculated intensities for the 
first magnetic peak (±1 ±1 1) are shifted to lower 20 
values in contradiction to the experiment. In contrast, 
the AFM1 structure perfectly describes the data if the 
magnetic moments are parallel to the a-direction. A fit 
with moments in b direction does not result in a correct 
description of the intensities. Even the free refinement 
of both a- and 6-components yields only negligible val- 
ues for the moment in 6-direction. This observation of 
the ordered moment along the a direction is in accor- 
dance with the measured magnetic susceptibility in Ref. 9 
which sharply drops along the a-direction below Tn ■ The 
AFM3 and AFM4 structures with c-components do not 
properly describe the magnetic intensities; see Fig. [31(a). 
In VOC1 the AFM1 magnetic structure with moments in 
a-direction is thus realized. In Ref. 9 , a mixture of the 
AFM1 and AFM2 was proposed basing on single-crystal 
diffraction, however, most likely, the twinning induced by 
the structural phase transition, which was not detected 
in that work, occurs in the used single crystal fully ex- 
plaining the observed intensities. 

The antiferromagnetic moment amounts to about 1.3 
fiB at 2 K which is in approximate agreement with the 
value reported in Refi^ but which is smaller than the full 
spin moment expected for the vanadium ion with a 3d 2 
electronic configuration. The value of the antiferromag- 
netic moment as a function of temperature is shown in 
Fig. [5] (h). The temperature dependency of this anti- 
ferromagnetic moment indicates a Neel-temperature of 
Tat -80 K. 
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FIG. 5: (color online) (a-d) The possible irreducible repre- 
sentations in space group P112/n for k = (1/2 1/2 1/2) with 
magnetic atoms in the primitive unit cell at (-0.25, 0.25, zo)- 
For (a) and (b) the presented spin direction along a indicates 
the direction of the ordered moment found in VOC1, which 
exhibits magnetic order according to the representation # 1. 



C. Magnetoelastic coupling in VOC1 

The structural refinements reveal a displacement of the 
vanadium ions away from the chlorine towards the oxy- 
gen sites in c-direction upon cooling, see Fig. [6] (a-d). 
The VO-planes become thus flatter at low temperature. 
In consequence the V-O-V bond angle along a increases 
upon cooling thereby enhancing the antiferromagnetic 
exchange interaction J a . The V-O-V path along a ex- 
hibits a bond angle considerably smaller than 180° but 
larger than the value of 90° where one expects the anti- 
ferromagnetic superexchange to vanisbZ ' 15 i 16 . In conse- 
quence the magnetic interaction J a depends sensitively 
on the bond angle and increases with its enhancement. 
The direct overlap of the V orbitals along b seem to dom- 
inate the magnetic interaction J\> which increases upon 
cooling due to the shrinking of the 6-parameter. The 
structural changes implied through the cooling appear 
thus to be driven by the magnetism which favors an en- 
hancement of the magnetic interaction parameters at low 
temperature. 

Magnetoelastic coupling is also evident in the coupling 
of the magnetic transition with the monoclinic distor- 
tion. In Fig. [6] (h) we plot the ordered magnetic mo- 
ment together with the monoclinic distortion as a func- 
tion of temperature demonstrating the close coupling. 
The structural phase transition from Pmmn to P112/n 
symmetry is of proper ferroelastic character, so that the 
monoclinic distortion is the primary structural order pa- 
rameter. 

More quantitatively, the preference of the AFM1 struc- 
ture can be understood with regard to the monoclinic dis- 
tortion of the structure. As noted above, the magnetic 
arrangement is antiferromagnetic along the V-V bonds, 
Jdiag (g reen in Fig. [2]), but ferromagnetic along Jf lag 
(red). Unfortunately, the neutron measurements suffer 
from the small vanadium scattering length causing en- 
larged error bars for variables which depend on one or 
on two V-position parameters. This effect is most visible 
in Fig. [6] (f) where the difference of the bond lengths of 
the diagonal V-V bonds (denoted by the red and green 
lines in Fig. ^ is shown. Therefore, an extremely long 
measurement has been performed on the D1A diffrac- 
tomcter in order to increase the statistics, and another 
synchrotron radiation experiment at 20 K has been per- 
formed at beamline B2 at Hasylab/DESY. In this syn- 
chrotron experiment several image-plate data collections 
with a total measuring time of more than one hour have 
been summed up. The results of these precise experi- 
ments are plotted together with the other neutron data 
in Fig. [6] (a-g). The error bars of the synchrotron mea- 
surements are of similar dimensions as the symbol size. 
Especially, the value of the difference of the diagonal V- 
V distances shown in Fig. [5] (f) has much smaller error 
bars and corroborates the neutron results. The structural 
parameters of this synchrotron measurement and of the 
neutron measurements at 2 K and at 150 K are listed in 
Tab. fl] The difference in the squares of the diagonal V-V 



5 



2 K: 


neutron 


D1A (ILL) 






SG 


a (A) 


b (A) 


c (A) 


T 


P112/n 


3.76755(10) 


3.29045(9) 


7.89582(20) 


90.21(0) 


Atom 


X 


y 


z 




V 


-0.25 


0.25 


0.11678(76) 




CI 


-0.25 


0.75 


0.33000(4) 




o 


-0.25 


0.75 


-0.04800(6) 






10 K: 


neutron 


D1A (ILL) 






SG 


a (A) 


b (A) 


c (A) 


7 


PU2/n 


3.76776(6) 


3.29076(5) 


7.89611(13) 


90.202(1) 


Atom 


X 


y 


z 




V 


-0.25 


0.25 


1 1 478('59 > l 




CI 


-0.25 


0.75 


0.33035(3) 




o 


-0.25 


0.75 


-0.04765(3) 




150 K: 


neutron 


D1A (ILL) 






SG 


a (A) 


b (A) 


c (A) 


T 


PmmnS 


3.77001(10) 


3.29656(10) 


7.90713(25) 


90 


Atom 


X 


y 


z 




V 


-0.25 


0.25 


0.12804(111) 




CI 


-0.25 


0.75 


0.32914(5) 




o 


-0.25 


0.75 


-0.04815(8) 




20 K: 


synchrotron 


B2 (DESY) 






SG 


a (A) 


b (A) 


c (A) 


7 


P112/n 


3.7654(3) 


3.2880(2) 


7.8894(5) 


90.207(1) 


Atom 


X 


y 


z 




V 


-0.25 


0.25 


0.11555(15) 




CI 


-0.25 


0.75 


0.33068(19) 




o 


-0.25 


0.75 


-0.04906(51) 





TABLE I: Results of the structural refinement of the neutron 
(D1A) and synchrotron (B2) measurements. 



distances is determined with much higher precision as the 
V z-parameter does not enter its calculation which only 
depends on the projection of the distances onto the a, b- 
planes, see Fig. |6](g). Since the x and y coordinates are 
fixed, the difference in the squares of the V-V distances 
is directly calculated from the lattice parameters. 

Indeed, the two diagonal V-V bonds (red/ green in Fig. 
O differ by 0.0073(15) A as can be seen in Fig. U (f). 
The splitting of the diagonal V-V distance is just a con- 
sequence of the monoclinic angle and clearly shows up in 
y/{d(V - VUag,!} 2 - {d(V - V) diaga Y, see Fig. ©(g). 
Compared to the orthorhombic structure where the diag- 
onal V-V bonds (red/ green in Fig. [5]) are equivalent, the 
monoclinic distortion renders the two antiferromagnetic 
exchange interactions inequivalent. We unambiguously 
find that the diagonal V-V distance associated with J^ iag 
is smaller than that associated with J diaQ which indicates 
a stronger orbital overlap accompanied by a larger ex- 
change interaction J^ iag thereby explaining the observed 
magnetic structure. Since the magnetic energy gain due 
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FIG. 6: Results of powder neutron diffraction (black circles) 
and synchrotron measurements (gray triangles), (a) V-O dis- 
tances for bonds in a- and 6-direction. (b) V-Cl distance, 
(c) V-V distances for V-ions within the same V-bilayer in 
diagonal direction (red and green bonds in Fig. [2}. (d) V- 
O-V bond angles for V-ions (with the same value of z) in 
a-direction (V-0-V) a and in 6-direction (V-0-V) 6 . (e) Cl-V- 
Cl angle and V-O-V angle in diagonal direction, (f) Ad(V- 
y)dia g denotes the difference of the two diagonal V-V dis- 
tances, (g) s/{d(V - V) dlag ,i} 2 - {d(V - V) dla3 ,2} 2 plotted 
with the monoclinic angle 7 (gray), (h) The antiferromagnetic 
moment (black) together with 7 (gray). (Lines are guide to 
the eyes.) 



to the lifting of the frustration is linear in the struc- 
tural distortion, and since the lattice energy cost scales 
quadratically with the distortion, one should always ex- 
pect such an effect if the magnetic coupling between two 
sub-units remains frustrated. However in most cases, the 
magnetoelastic coupling would be too weak to generate 
an observable effect. The case of VOC1 seems to profit 
from the fact that the diagonal interaction is modulated 
via the V-V distances without variation of the V-0 dis- 
tances. The structural modulation in VOC1 causes thus 
less lattice energy loss than for instance a direct bond- 
distance modulation. 

The splitting of the diagonal V-V distances is accom- 
panied by a small splitting of the V-O-V angles leading 
to a smaller angle at the antiferromagnetic bond which 
seems to contradict the first estimation that the magnetic 
interaction should increase with the bond. However, in 
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this tig system the direct overlap varied through the di- 
agonal V-V distance seems to be the dominant effect. 

IV. CONCLUSION 

In conclusion we have studied the nuclear and mag- 
netic structure of VOC1. Synchrotron diffraction mea- 
surements reveal a structural phase transition Pmmn — > 
P2/n which accompanies the magnetic ordering in VOC1. 
By neutron diffraction we were able to refine the mag- 
netic structure within a monoclinic 2x2x2 unit cell and 
derived a 3D-antifcrromagnetic structure with an antifer- 
romagnetic moment pointing in a-direction. At 2 K this 
moment amounts to about 1.3 hb which is a significantly 
reduced value for a V 3+ -ion in this 3d 2 -system. 

Magnetoelastic coupling plays an important role in 
VOC1 due to the fact that the magnetic exchange inter- 
action depends sensitively on the V-O-V bond angles and 



on the orbital overlap. The modulation of the bond an- 
gles and of V-V distances involves less lattice-energy cost 
than e.g. bond-length modulation, so that rather strong 
structural effects occur. First the VO layers flatten upon 
cooling which can be well understood in terms of an en- 
hanced magnetic interaction at low temperatures. The 
anisotropic thermal expansion further strengthens the en- 
hanced magnetic interaction at low temperature. Second, 
the monoclinic distortion accompanying the magnetic or- 
dering lifts the frustration of the coupling between two 
magnetic subsets as the diagonal V-V distances shrink in 
one diagonal direction enhancing the orbital overlap and, 
thus, the exchange interaction. 
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